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ABSTRACT 

The spin is an important but poorly constrained parameter for describing supermassive black holes (SMBHs). 
Using the continuity equation of SMBH number density, we explicitly obtain the mass-dependent cosmological 
evolution of the radiative efficiency for accretion, which serves as a proxy for SMBH spin. Our calculations 
make use of the SMBH mass function of active and inactive galaxies (derived in the first paper of this series), the 
bolometric luminosity function of active galactic nuclei (AGNs), corrected for the contribution from Compton- 
thick sources, and the observed Eddington ratio distribution. We find that the radiative efficiency generally 
increases with increasing black hole mass at high redshifts (z > 1), roughly as 77 cx M^^, while the trend reverses 
at lower redshifts, such that the highest efficiencies are attained by the lowest mass black holes. Black holes 
with M, > IO^'^Mq maintain radiative efficiencies as high as « 0.3-0.4 at high redshifts, near the maximum 
for rapidly spinning systems, but their efficiencies drop dramatically (by an order of magnitude) by z w 0. 
The pattern for lower mass holes is somewhat more complicated but qualitatively similar Assuming that the 
standard accretion disk model applies, we suggest that the accretion history of SMBHs and their accompanying 
spins evolve in two distinct regimes: an early phase of prolonged accretion, plausibly driven by major mergers, 
during which the black hole spins up, then switching to a period of random, episodic accretion, governed by 
minor mergers and internal secular processes, during which the hole spins down. The transition epoch depends 
on mass, mirroring other evidence for "cosmic downsizing" in the AGN population; it occurs at z « 2 for 
high-mass black holes, and somewhat later, at z w 1, for lower-mass systems. 

Subject headings: black hole physics — galaxies: evolution — quasars: general 



1. INTRODUCTION 



It h as been realized since the pioneering work of ISoltarJ 
( [19821) that supermassive black holes (SMBHs) located at 
the centers of galaxies assemble their mass predominantly 
through accretion, an inference that has been further rein- 
forced fro m consideration s of the cosmic X- r ay background 
jFabian & Iwasawal 119991: lElvis et al.1 120021: iMarconi et all 
I2OO4.) . However, how SMBHs are fueled remains an outstand- 
ing unsolved issue. 

The spin of SMBHs traces the angular momentum of the 
accreted material. As such, it can serve as a powerful cos- 
mic probe of SMBH feeding. However, spin is highly elusive 
to measure since its general relativistic effects emerge in the 
very vicinity of the horizon (typically within a few tens grav- 
itational radii). According to standard accretion disk theory, 
the radiative efficiency of energy conversion is closely linked 
to black hole spin through the marginally stable orbit, which 
is a function of the spin. The binding energy of the material 
in the accretion disk is locally radiated, and, after crossing the 
marginally stable orbit, the material freely falls into the hole 
without losing further energy due to the torque-free condition 
there. As a result, the total amount of energy converted into 
radiation is the bin ding energy b etween the marginally stable 
orbit and infinity ( lThornelll974l) . Specifically, the radiative 
efficiency incre ases monotonically with black hole spi n (e.g., 
see Figure 1 of iMartmez-Sansigre & Rawlingsll201 lal) . This 
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link allows us to analyze the net angular momentum of the ac- 
creted gas by quantifying the radiative efficiency, and hence 
obtain clues on ho w SMBHs are fed. By presuming a fixed 
radiative efficiencv. lSoltanI (119821) connected the mass growth 
rate of SMBHs with the luminosity function (LF) of active 
galactic nuclei (AGNs). Subsequent studies explored the cos- 
mic growth of SMBHs by comparing the local SMBH mass 
density w ith the integrated energy density of AGNs across 
time (e.g.,|C hokshi & Turner" 1992'; 'Small & Blandford"1992|; 
lYu (& Treinaine 2002; M arconi et al.. ,2004; Shankar et all 
2004; Cao & Li 2008; Ca(^(2010|; see also a review of |S|iankai| 
l2009i) . These studies found that an average radiative ef- 
ficiency of r] K, Q.\ yields a local mass density consistent 
with observational constraints. From the theoretical point of 
view, if the angular momentum of the accreted material stays 
aligned with the spin axis of of the black hole and the direc- 
tion remains unchanged, the black hole will be rapidly spun 
up to a maximally ro tating system and attain a radiative ef- 
ficiency of 1] w 0.42 (lThorndll974l) . The inconsistency be- 
tween this estimate of 77 and that based on Soltan's argument 
potentially intimates that episodic transitions of angular mo- 
mentum of accreted material during the active phases of black 
holes (i.e. random accretion) may play an important role in 
the cosmological evolution of SMBHs (e.g.. lKing et al.ll2008t 
IWang et al.l2009h . It is thus expected that the spin of SMBHs, 
and hence their radiative efficiency, evolves with redshift. 

,Wang et al., (.2009i) constructed a formalism to determine 
the evolution of the radiative efficiency, which depends com- 
pletely on observables. When applying this formalism to sur- 
vey data, they found that SMBHs are spinning down with cos- 
mic time since z < 2. Such an evolutionary trend of spins 
suggests that SMBH growth is driven by episodic random ac- 
cretion. Since different SMBH populations may undergo dif- 
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ferent accretion histories, it is important to verify whether the 
spin evolution depends on black hole mass. This can give fur- 
ther insights into how SMBH activity is triggered and how the 
angular momentum of the accreted material ultimately influ- 
ences SMBH spin. 

Deep surveys in the X-rays and in other bands have es- 
tablished that AGN activity exhibits cos r nic downsizing (e.g. , 
Ueda et al.ll2003t iBongiorno et al.ll2007HHopkins et al.ll2007l: 



Cirasuolo et al. I|2010!). GeneraUy, the space density of AGNs 



with low luminosity peaks at lower redshift than that of AGNs 
with high luminosity. Two distinct scenarios can result in such 
luminosity evolution. If black holes are accreting at near- 
Eddington luminosity once they become active, downsizing 
can be caused by activity s hifting toward low-mass black 
holes at low redshift (e.g . . iHeckman et al] l2004t IShankad 
I2009t ISchulze & Wisotzka |201(]|) . Alternatively, downsiz- 
ing can result from a decrease o f the average accr etion rate 
onto black holes at low redshift (iBabic et al.ll2007h : in other 
words, the entire black hole population shines with high Ed- 
dington ratio at high redshift and then slowly fades out over 
time. These two scenarios can be tested using information 
on black hole demographics. Studies of local AGN samples 
unambiguously find that local active black holes are typically 
an order of magnitude less massive th an the typical inactive 
black holes residing in normal galaxies ( IHeckman et al.ll2004l : 
iGreene & Holl2007l) . Using the virial method to measure the 
black hole m asses of SDSS quasars out to redshift z w 4, 
iLabita et"an (2009a b) show that the maximum mass of the 
active black hole population notably increases with redshift. 
These observations seem to suggest that cosmic downsizing 
arises from activity shifting over mass. Recalling that the an- 
gular momentum carried into the black hole along with mass 
accretion changes the black hole spin, we may expect that the 
spin evolution of black holes would follow a behavior similar 
to that of the redshift evolution of AGN activity. 

Assuming that SMBHs grow mainly through accretion, we 
extend the study of Wang et al. (2009) to explicitly show the 
mass-dependent behavior of the cosmic evolution of SMBH 
spin. In Section 2, we construct a generalized equation that 
expresses the radiative efficiency as a function of black hole 
mass and redshift. Section 3 shows the Eddington ratio dis- 
tribution used in our calculations, and Section 4 derives the 
SMBH mass function of AGNs, including Compton-thick 
sources. We then present the results for SMBH growth in 
Section 5 and SMBH spin evolution in Section 6. The impli- 
cations of our results on accretion scenarios are discussed in 
Section 7. Conclusions are summarized in Section 8. 

Throughout the paper, we adopt a cosmological model with 
nM = 0.3, 17a = 0.7, and //o = 70km s-1 Mpc"'. 

2. THE CONTINUITY EQUATION OF BLACK HOLES 

Let N(t,M,) be the SMBH mass function, including both 
active and inactive black holes, which specifies the number 
of black holes per unit comoving volume and per unit mass 
at cosmic tim e t. Its evolution is descr ibed by a continuity 
equation (e.g.. lSmall & Blandfordll 19921) 



dN(t,M,) d 



dt 



+ -^[N{t,M,){M,)]=S{t,M,), 



(1) 



where (M,) is the mean mass accretion rate for inactive 
and active SMBHs with mass M, at time r, and 5(r,M,) 
is the source term that accounts for black hole mergers. 
As in most previous works (e.g., .Small & Blandford. 19921 : 



lYu & Tremain3l2002t iMarconi et al.ll2004l). we neglect black 
hole mergers, setting SUM.) = 0. IShankar et all ( I2009ll2010l) 

investigated the importance of black hole mergers on the 
evolution of the SMBH mass function and concluded that 
the effect of mergers is minor comp ared with mass a ccre- 
tion. Numerical s imula tions by iVolonteri et alJ (l2005h and 
iBerti & Volonteril (120081) also verified that mass accretion 
dominates over mergers in determining the mass growth and 
spin distribution of black holes. 

The radiation power of a SMBH with mass accretion rate 
Mace is -L = jyMaccC^. Because of radiative losses, the mass 
growth rate of the SMBH is M, = (1 -?7)Macc, yielding M, = 
(1 -rfjXL^ii/rjc-. Here the Eddington ratio is defined as A = 
L/LEdd, where Ledd = 4-TTGM,mpc/<j'Y, G is the gravitational 
constant, nip is the proton mass, c is the light speed, and aj 
is the Thomson electron scattering cross section. With the 
help of the duty cycle of SMBHs, S{t,M,), the mean mass 
accretion rate can be written as 



{M,)=S{t,M.)M, = S{t,M,) 



1-7? A(f,M.)LEdd 

7] C 



(2) 



where X{t,M,) is the mean Eddington ratio determined from 
observations (see Equation [121 below) . The duty cycle is usu- 
ally defined as 



5(f,M.): 



Nc(t,M,)+NACN(t,M.) 



(3) 



where Nc(t,M,) and A^AGN(f ,^.) are the SMBH mass func- 
tions of galaxies and AGNs, respectively, and, accordingly 



N(t,M,) = NG{t,M,)+NAGNit,M,). 



(4) 



Combining the above equations, we rewrite the continuity 
equation 



dN(z,M») dt d 



dz 



dzdM, 



1 -77 A(z,M.)LEddA'AGN(z,M.) 



(5) 

where and hereinafter we substitute the cosmic time t with 
the corresponding redshift z- Integrating Equation (|5]) over 
M'^ fromM, to 00, we obtain the radiative efficiency 

7?-'(z,M.)= 1 + — — - - / N{zM',)dM',, (6) 
m(z,M,) \dzj ozJm. 



where the AGN luminosity density is 

u(z,M,) = A(z,M.)LEddA^AGN(z,M.), 



(7) 



and we use the boundary condition u(z,M, — > 00) = 0. Equa- 
tion (|6]l is the generalized 7]-equation of Wang et al. (2009). 
The underlying rationale of this generalized 77-equation lies in 
the assumption that mass accretion only increases the SMBH 
mass, but keeps the SMBH number density conserved. Sub- 
sequently, variations of the SMBH number density in the in- 
terval [M,,oo] should come from the "inflow" of SMBHs 
through the boundary = M, due to mass accretion. As 
a result, we obtain the radiative efficiency for SMBHs with 
any given mass of M, from the AGN luminosity density. 

To solve Equation (|6]l, it is imperative to obtain three in- 
gredients that appear on the right-hand side: the SMBH mass 
functions of normal galaxies and AGNs, and the E ddington 
ratio distribution. In the first paper of this series dLi et al.l 
I20T1I Paper I), we have derived the SMBH mass function of 
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Fig. 1 . — SMBH mass function of AGNs. Solid and dot-dashed lines are the calculated mass function using Equation (To) with a log-normal and power-law 
distrib ution of Eddington ratios, respectiv ely. The s haded areas represent a typical error of A log^AON ~ ±0 . 2 dex from the AGN bolometric LF of Hopkins et al] 
I2007D . Asterisks in the 2 = bin are from lGreene & Ho ( 2007). and squares are from[Vestergaard & Osmeil 120091) . who also compiled the SMBH mass function 
at 0.025 < ; < 0.5 from the Bright Qu asar Sur vey. For the sake of a comparison, we superpose it in the z = bin. Note that the turnover toward the low-mass end 
of the mass function of|^stergaard & Osmeil is due to high survey incompleteness jKelly et al. 2013). 



normal galaxies out to z w 2 using the latest galaxy luminos- 
ity and stellar mass functions. In what follows, we show in 
detail how to obtain the SMBH mass function of AGNs from 
the observed Eddington ratio distribution. 



thereby appearing as low-luminosity AGNs (lHoll2009allbl) . 

With the above lines of observations, we employ two types 
of Eddington ratio distribution in our calculations. We express 
the log-normal distribution as 



3. EDDINGTON RATIO DISTRIBUTIONS 

Observations thus far have not yet reached a consensus 
on the observed Eddington ratio distribution, especially at 
high redshifts. The challenges arise from systematic uncer- 
tainties in measuring black hole mass using the virial rela- 
tion and the flux limit of surveys. However, among ear- 
lier studies, two types of Eddington ratio distributions are 
preferred. Generally, the observed Eddington ratios de- 
rived from the brig ht AGN (quasar) samples exhibit a log- 



rived from the brig ht AGN (quasar) samp ! 
normal distribution (Kollmeier et al.ll2006t 
Kauffmann & Heckman 2009; Ke llv et al.l 



Shen etal.ll200i 



Eauffmann & Heckman 2009; Ke lly et al.l l2010h . Bv con 
trast, those derived from samples that include fainter sources 
or low-luminosity AGNs display a power-law distribution 
and/or an additional l og-normal componen t in the ^ high- 
Eddington ratio regime (Heckman et al."'2004'; 'Yu et al.' 2005; 
Hopkins & Hernauist 2009; Kauffmann & Heckman 2009; 
Schulze & Wisotzkill2010l) . In particular, the peak and dis- 



persion of the log-normal distribution are found to be almost 
independent of lu minosity and redshift (iKoUmeier et al.l2006l: 
IShen et aTll2008l) . Through analysis of the properties of the 
host g alaxies selected from SDSS, iKauffmann & HeckmanI 
(12009*) argued that in the log-normal regime black holes self- 
regulate their growth in environments with a plentiful supply 
of cold gas, so that the Eddington ratio is retained at a uni- 
versal rate independent of luminosity (black hole mass) and 
the properties of the host galaxies. By contrast, in the power- 
law regime the gas has run out and the central black holes are 
probably being fueled by mass loss from stars in the bulges, 



P{\)d\og\ = 



exp 



(logA-M)' 
2a'- 



c/logA, (8) 



where /i is the logarithm of the peak Eddington ratio and a is 
the dispersion. Below we will show how to determine fi. The 
dispersion a is found to be inessen tial to our results, and we 
set a = 0.3 (iKoUmeier et al.ll2006h . The power-law distribu- 
tion is given by 



P(A)<ilogA = Co 



exp 



"Ao 



c/logA 



A > An 



(9) 

where Co is the normalization, Ao is the characteristic ra- 
tio, K is the power-law index, and Am m is the lo wer cut-off 
(iHopkins & H ernquist 200§; see also Cad 120 1 0|) . Accord 



Hopkins & HernquistI 



ing to their feedback-regulated model, 
(120091) suggest typical values of Ao « 0.2-0.4 and k « 0.3- 
0.8. In our calculations, we set Ao = 0.3 and k = 0.6 as fiducial 
values. Note that for the power-law distribution, a lower cut- 
off Amin is required to completely determine the normalization 
Cq. We adjust the values of Amin so that the calculated SMBH 
m ass function of A GNs matches the observed one. 

ICao & Lil (I2008h derived the mean Eddington ratio for a 
given black hole mass and redshift from the observed Edding- 
ton ratio distribution. For clarity, we present the details here. 
The SMBH mass function of AGNs is calculated by combin- 
ing the observed Eddington ratio distribution and the AGN 
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Fig. 2. — Mean Eddington ratio as a function of (top) black hole mass and 
(bottom) redsliift. 

bolometric LF as 

A^AGNfcM.)= / $(z,Lboi) '^'°^^'°' f(AVlogA, (10) 
] dlogM, 

where $(z,Lboi) is the AGN bolometric LF. Applying Bayes' 
theorem, the Eddington ratio distribution for a given black 
hole mass M, is 



Fig. 3. — Duty cycle as a function of (top) black hole mass and (bottom) 
redshift. 



a;(z,A|M.) = 



(11) 



A^agn(z,M.) 

The mean Eddington ratio of AGNs with M, at redshift z is 

Xz,M,) = j \uj{z,X\M,)d\og\ (12) 

where / w(z, A|M,)(ilogA = 1. 

In Figure [T] we compare with observations the SMBH 
mass function of AGNs calculated by Equation ( fTOl i using 
two types of Eddington distributions. To delineate the evo- 
lution of the AGN populations, we adopt Hopkins et al.'s 
(2007) bolometric LF, which combines a large set of AGN LF 
measurements in the rest-frame optical, soft and hard X-ray, 
and near-IR and mid-IR bands, spanning a range of bolomet- 
ric lum inosities from ^ 10"*^ to lO'*^ ergs"'. iGreene & Hoi 
( 12007b measured the local SMBH mass function of broad- 
line AGNs from SDSS using standard virial relations. Based 



on the same method, IVestergaard & Osmed (l2009h compiled 
the SMBH mass function of broad-line AGN samples drawn 
from the Large Bright Quasar Survey, the Bright Quasar Sur- 
vey, and SDSS, covering a range of redshift from the local 
epoch up to z = 5. Note that t he turnover toward the low-mass 
end of the mass function of IVestergaard & Osmei* is due to 
the incompleteness of the surveys. iKell v et al. (2010) care- 
fully estimated the incompleteness of the SDSS sample of 
5festergaard & Osmer and found that it is highly incomplete 
for M, < 10'* M© at z > 1. Since these measurements just 
cover the unobscured population of AGNs, a correction factor 
has to be applied to Equation (fTOl i to account for such a se- 
lection bias before performing this comparison (see Section 4 
below). By fine tuning the value of ^ for the log-normal dis- 
tribution of Eddington ratio, we find that = -0.6 produces an 
AGN SMBH mass function that is in excellent agreement with 



the ob served ones, particularly for the data of iGreene & Hoi 
(120071) at z = 0, as shown in Figure [T] Similarly, a value of 
Amin = 0.1 for the power-law distribution also gives an almost 
identical mass function. However, it is probably unphysical 
that the Eddington ratio for the power-law distribution is cut 
off at Amin = 0.1. Nevertheless, whichever type of Edding- 
ton ratio distribution is used has very little influence on the 
final results in the sense that the resultant AGN SMBH mass 
functions are very insensitive to the actual choice (see Figure 
[T]). Hereafter, we only use the log-normal distribution in our 
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Fig. 4. — Luminosity density li (see Equation|7) as a function of redshift, 
for different black liole masses. 

calculations. 

Figure|2]shows the dependence of the mean Eddington ratio 
on black hole mass and redshift. We find that the mean Ed- 
dington ratio roughly spans values in the range ~ 0.1 -0.3, 
slightly increasing with redshift and decreasing with black 
hole mass. These trends between the mean Eddington ratio 
and black hole mass and redshift are due to the power-law 
shape of the AGN LF and the redshift evolution of the slope 
of the power law. At the same time, we know that nearby 
massive active galaxies generally contain black holes accret- 
ing at an Eddington ratio substantially below unity (e.g., iHol 
|2p08, 2009a). This does not conflict with the present results, 
considering that local luminous AGNs are quite rare. More- 
over, we here focus on efficient accretion phases with high 
Eddington ratios, during which black hole growth predomi- 
nately happens (e.g., Hopkins et al. 2006; Xu & Cao 2 010,) . 

4. ACCOUNTING FOR COMPTON-THICK AGNS 

Given the Eddington ratio distribution and AGN bolomet- 
ric LF, it is trivial to calculate the SMBH mass function 
for AGNs using Equation ( fTOb . However, it is important to 
note that different bands are subject to diffe rent selection bi- 
ases. Optical LFs miss obscured AGNs (e.g.. lBongiorno et al.l 
l2007h . while hai'd X -ray surveys (~ 2- 10 keV) can trace only 
AGNs that are not completely Compton-thick (e.g.. lHasingeii 
2008h. namely sources with column densities A^h ^ 10 cm 



Treister et alj ( 120101) recently presented an evolution model 



for AGNs hosted by galaxies undergoing major mergers and 
find that Compton-thick sources play an important role in the 
mass growth of SMBHs. Meanwhile, evidence is mounting 
from deep X-ray surveys, in combination with IR surveys, that 
the population of Compton-thick AGNs with intermediate lu- 
minosities is of the same order as that of Compton-thin AGNs 
(sources wifli A^h » 10^^- 10^"^ crn"^; e g..lDaddi et al.ll2007l: 
[Alexander etaTIIIOOSl: iFiore et al.ll2009HTreister et al.ll2009t) . 
Population synthesis models of the cosmic X-ray background 
also predict the existence of a large number of Compton-thick 
AGNs (e.g., Gilli et al. 2007; Draper & Ballantyne 200^)- 

How the fraction of Compton-thick AGNs evolves with 
redshift is not well known. Even the overall evolution of 
Compton-thin AGNs is in debate. After correcting for se- 
lection bias in their hard X-ray-sele cted AGN sample s pan- 
ning the redshift range z = 0-4. iTreister & Urrvl (|2006|) find 



Fig. 5. — Evolution of black hole mass (dot-dashed lines) and associated 
radiative efficiency (solid lines) as a function of redshift. The initial masses 
are set to 10^ lO' ^ 10^ and 10^^ Mq at redshift zo = 2.5. The symbols, 
squares and circles, indicate when the black hole mass has e-folded once and 
twice, respectively. 

that the fraction of obscur ed AGNs increase s with redshift as 
(l+zT, with a w 0.3-0.5. iBallantyne et alldJoOei) . modeling 
the X-ray background, also report a similar evolution trend of 
(l+z)"^. However, other studies cast do ubt on the redshift 
evolution of th e obsc ured AGN fraction (|Ueda et al. 120031: 
iDwellv & P age 2006; GilH et al. 2007; Lama stra et al.lboOSh . 
Based on their AGN population model, Gilli et al.l ( 120101) ar- 
gue that by including a proper /T-correction an intrinsically 
constant obscured AGN fraction can also represent the data. 

Motivated by the above results, we parameterize the red- 
shift dependence of the fraction of Compton-thick for all 
AGNs as 

/c = /c,o(l+zr, (13) 

with fcfi ~ 0.3 and a « 0.3. Here we determine /c,o 
by assuming that the local fraction of obscured (Compton- 
thin and Compton-thi ck) to unobscured AGNs is ~ 3 : 1 
( ITreister & Urrvll2006l) . and that there is an equal abundance 
of Compton-thin and Compton-thick sources. As described 
later, our final results are quite insensitive to the value of 
a; the refore, we set a = 0.3, consistent with ITreister & Urrvl 
( 120061) and Ba llantvne et al.l (120061) . On the other hand, it 
is well established that the obscured frac tion of Compton- 
thin AGNs decrease s with luminos ity (lUeda et al.l 120031: 
iLa France et al.l [20051: iHasingeJ 120081) . The dependence of 
obscuration on luminosity for Compton-thick AGNs, how- 
ever, is still poorly understood. We neglect this effect, but 
we discuss the implications in Section 6.2.3. 

5. EPISODES OF SMBH GROWTH 

With the AGN SMBH mass function determined from the 
previous sections, we can calculate the AGN duty cycle, as 
shown in Figure |3] At low redshifts, the duty cycle decreases 
with increasing black hole mass, indicating that the fraction of 
galaxies that are active also declines with mass (iGreene & Hoi 
12007; Shankar 2009). Conversely, as a consequence of the 
cosmic downsizing of AGN activity, at high redshifts the duty 
cycle reverses trend and increases with black hole mass. Duty 
cycles for lO^M© and lO^M© black hole are 5 w 4 x lO"-' 
and 8 X lO"'' at z = 0, respectively, but rise up to ^ 0.06 and 
~ 0.5 by z = 2. As shown in the bottom panel of Figure[3] the 
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Fig. 6. — Radiative efficiency as a function of black hole mass at different redshifts. Solid lines and dot-dashed lines are the efficiency using the SMBH mass 
function derived from the galaxy luminosity and stellar mass function, respectively. Shaded areas denote typical en'ors of A log = ±0.5 dex from SMBH mass 
function and AGN LF. 



active fraction of all black holes increases significantly from 
~ 10"^ at z = to 0.1 at z = 2, in concert with the feature 
of AGN LF that displays a rapid rise up till z = 2 (see also 
iWanpet al.ll20()6h . 

With the duty cycle determined, the active time of AGNs is 
simply given by t(z,M,) = (5(z, M,W(z)~', whe re H{z) is the 
Hubble parameter at redshift z (lPeacocklll999h . As a result, 
the active time of AGNs is of the order of ~ 10^ yr in the 
local Universe and ~ 10** yr at z = 2. This is consist ent with 
the ob servational estimates from other methods (see iMartinil 
|200jfor a review). 

Figure |4] illustrates the AGN luminosity density u(z,M,), 
defined by Equation (jTji, and its variation with redshift for 
different black hole masses. Since the mean Eddington ratio 
is approximately constant over black hole mass and redshift, 
u(z,M,) just reflects the evolution of the AGN SMBH mass 
function. Again, the behavior of cosmic downsizing is evi- 
dent; more massive active black holes reach their maximum 
number density earlier than less massive ones. 

With the mean accretion rate (M,), we can trace the SMBH 
growth history since redshift zo as 



quantified by the Salpeter (1964) time 



1] 



{l-ri)X4TTGmp 



= 2.3 X 10** 



lOry 



1- 



A 

02 



yr. 

(15) 

It is obvious that the e-folding time scales in proportion to 
the radiative efficiency. The decline of radiative efficiency 
makes black hole growth possible, especially for low-mass 
black holes at low redshift, considering that the AGN life time 
is of the order of 10^ yr at that epoch (see also the discus- 
sions of iKing & Pringlell2006t iNetzer & TrakhtenbrotI 120071: 
iNetzer et al.ll2007b . 

An intriguing issue that has attracted much attention over 
the years is what is the mass range of seed black hol es and 
how they g row to prese nt-day SMBHs (e.g., Volont eri et alJ 
|2008; Wan g et aLll2008h . In principle, an extrapolation of the 
present approach to much higher redshifts, say z > 10, may 
be helpful to place constraints on the properties of seed black 
holes. This is beyond the scope of the present paper, and we 
defer this investigation to future work. 



6. SMBH SPINS 



M,(z)=Mo + 



{M,[z',M.(z')])-^dz', 



(14) 



6.1. 



where Mq is the initial black hole mass at redshift zo- Figure|5] 
plots the growth history of black holes with initial masses 10^, 
10^ ^ 10^ and lO**-^ Mq at zo = 2.5. The symbols (squares 
and circles) indicate when the black hole mass has e-folded 
once and twice. For a hole with an initial mass of IO^Mq, 
it grows to ~' 1.5 X 10^ Mq at z ~ 0.3, e-folding roughly 3 
times. A 10^-^ Mq hole grows just to ^ 1.0 x 10^ at z ~ 
0.3, e-folding only once. Recall that the e-folding time is 

* Our estimates of e-folding times and mass growth rates are quahtatively 
consi stent with previous studies (e.g., Marconi et al. 2004; Merloni & Hein3 
12003) . Note, however, that these studies assumed a constant radiative effi- 



Evoliition of the Radiative Efficiency 

In Figure |6] we present the radiative efficiency, calculated 
using Equation as a function of SMBH mass, using the 
SMBH mass function derived from both the galaxy LF and the 
galaxy stellar mass function. Shaded areas represent typical 
errors of A log 77 = ±0.5 dex from the SMBH mass function 
(- 0.3 dex) and AGN LF (- 0.2 dex). We find that at low red- 
shifts (e.g., z ~ 0.3) SMBH of all masses are accreting mate- 
rial with a relatively low radiative efficiency. At high redshifts 
(e.g., z ~ 1-2), although the uncertainties are large, the ra- 
diative efficiency generally increases with black hol e mass. A 
rough correlation is 77 cx M2, with 7 w 0.5. Indeed. ICaom 
(|2008i) argued on the basis of the Soltan argument that the 

ciency. 
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Fig. 7. — Radiative efficiency evolution for different black hole masses. 
Solid lines and data points are the efficiency using the SMBH mass function 
derived from the galaxy luminosity and stellar mass function, respectively. 
Shaded ai'eas and eiTor bars denote typical errors of A log r] = ±0.5 dex from 
the SMBH mass function and AGN LF. Vertical dotted lines represent the 
redshift below which the efficiency begins to decline. 

radiative efficiency should increase with black hole mass so 
as to guarantee that the mass function of local relic SMBHs 
matches the mass function of galaxies. Previous studies, as- 
suming a fixed radiative efficiency ov er black hole mass and 
redshift, co ncluded that ?/ w 0.1 (e.g.. lYu & Tremain3l2002t 
iMarconi eTa l. 2004). The values of the radiative efficiency 
we obtain here are consistent with this average value. 

The redshift dependence of the radiative efficiency is plot- 
ted in Figure |7] for different black hole masses. Solid lines 
show the results obtained using the SMBH mass function de- 
rived from the galaxy LF, while the superposed data points 



make use of the SMBH mass function derived from the galaxy 
stellar mass function, computed for specific redshift bins. 
Again, shaded areas and error bars represent typical errors 
of A log 77 = ±0.5 dex. For high-mass black holes (M, > 
10*^^ M0), the radiative efficiency maintains a high value of 
77 w 0.3 at z w 2, but then strongly decUnes toward lower red- 
shifts. The radiative efficiency for M, « 10^ Mq at z w 2 
marginally exceeds the maximum efficiency of w 0.4 al- 
lowed for extreme Kerr black holes. This may arise from 
the unrealistic treatment of the mean Eddington ratio distri- 
bution beyond a black hole mass of IO'^Mq. Another cause 
may be an overestimate of the correction factor for Compton- 
thick AGNs in Equation ( fT3b , if the factor is generally anti- 
correlated with luminosity. For low-mass black holes, the 
evolution is more complicated. The radiative efficiency ap- 
pears to plateau at z w 1-2 and th e n dec reases dramatically 
toward z = 0. Recentlv. lWang et"al] ( 120091) obtained the aver- 
age radiative efficiency for black hole masses 
with redshift z < 2. They found a similarly rapid decrease of 
?7 from z ~ 2 to the local epoch (see their Figure 2a). 

Figure|5]shows how the black hole mass and the associated 
radiative efficiency evolve over time. Interestingly, there ex- 
ists a peak in the efficiency curve: for a given initial mass 
at redshift z = 2.5, the efficiency gently rises to a local max- 
imum, and then sharply falls off. The rise of the efficiency 
during black hole growth is due to the positive correlation be- 
tween efficiency and mass, as shown in Figure |6l while the 
subsequent decline is due to the systematic decrease of the ef- 
ficiency toward low redshifts. It seems that the efficiency evo- 
lution can be characterized by two regimes, one during which 
the efficiency increases and another during which it decreases. 

6.2. Evolution of SMBH Spins 

If the standard accretion disk model applies to the AGNs 
under consideration, we can elucidate their spin evolution 
through the radiative efficiency obtained above (e.g.. lThornel 
[19741) . Since black holes gai n their mass predom i nately dur- 
ing th eir quasar phase (e.g., iHopkins et alj|20()6t IXu & Caol 
120 101) . the standard accretion disk model is quite a reasonable 
approximation. From Figure |6l it is apparent that at high red- 
shifts (z > 1) black hole spin increases with mass. At low 
redshifts the dependence of spin on mass seems reversed, al- 
though large uncertainties prohibit us from reaching a firm 
conclusion. Volonteri et al. (2007) have drawn a connection 
between black hole spin and galactic morphology and argued 
that SMBHs in elliptical galaxies possess higher spins than 
those in spiral galaxies. Such a morphology-related spin dis- 
tribution has been used to explain the observed radio- loudness 
bimodality in nearby AGNs (e.g.. lSikora et al.|[2007b . Our re- 
sults at z ~ 0.3-0.6 are not qualitatively consistent with this 
scenario. One possible explanation for this inconsistency is 
that radio-loud galaxies comprise only a tiny fraction of the 
galaxy population, whereas our results pertain to a global av- 
erage of all galaxies. Moreover, our current approach is un- 
able to constrain the dispersion of the spin distribution at any 
given redshift or mass. 

Most importantly. Figure |7] shows that high-mass black 
holes begin to spin down earlier than low-mass black holes. 
To guide the eye, the vertical dotted lines in the figure mark 
the redshift below which the efficiency begins to decline. Re- 
calling the cosmic downsizing of AGN activity, there appears 
to be a connection between black hole spin and accretion. In- 
deed, if random accretion plays a role in black hole growth. 
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Fig. 8. — Influences of the observed Eddington ratio distribution on radia- 
tive efficiency for black holes with mass M. = 10**^ Mq . The dependence on 
the mean (p) and dispersion (a) of the distribution is illustrated in the top and 
bottom panels, respectively. 

spin evolution inevitably follows cosmic downsizing just as 
AGN activity does, since high-mass black holes become ac- 
tive earlier and hence are spun down earlier by random accre- 
tion. However, we observe two regimes in the evolution of 
the radiative efficiency: SMBHs first spin up, and then down. 
This indicates that more complex accretion scenarios need to 
be explored. We return to this point in Section 7. 

6.3. Influence of Uncertainties 

We explore the reliability of the radiative efficiencies in 
light of some of the assumptions used in our calculations. 

6.3.1. SMBH Mass Function of Galaxies 

The SMBH mass function is a crucial input because it deter- 
mines the mass growth rate in Equation (|6]l. In Paper I, we de- 
rived the SMBH mass function of galaxies out to redshift z ~ 2 
using the latest luminosity and stellar mass functions of field 
galaxies to constrain the masses of their spheroids. SMBH 
masses were inferred through a locally calibrated empirical 
correlation between black hole mass and spheroid mass. The 
SMBH mass functions derived from the galaxy luminosity 
and stellar mass functions show very good agreement, both 
in shape and in normalization, testifying to the robustness of 
our results. After carefully examining various sources of un- 
certainties, we found that the total uncertainties on the SMBH 
mass function are generally within ^ 0.3 dex, and possibly 
even smaller for lower mass black holes (see Paper I for de- 
tails). As a consequence, uncertainties in the SMBH mass 



Fig. 9. — Influences of the fraction of Compton-thick AGNs on radiative 
efficiency for black holes with mass M. = 10*^ Mq . The dependence on the 
nomialization (fc.o) ™d redshift dependence (a) of the Compton-thick frac- 
tion (see Equation 13) is illustrated in the top and bottom panels, respectively. 



function introduce an uncertainty less than 
resultant radiative efficiency. 



0.3 dex to the 



6.3.2. Eddington Ratio Distributions 

The observed distribution of Eddington ratios is described 
as lo g-normal, independe nt of lu minosity and redshift (e.g., 
Koll meier et al.i r2006: Sh en et ani200 8). Such a distribution 
may be realistic in the context of a self-regulation model 
for accre tion of black holes embedded in a gas-rich envi- 
ronment dKauffmann & HeckmanI [20091) . To assess its im- 
pact on our results, we calculate the radiative efficiency for 
three values of the mean Eddington ratio (logarithmic val- 
ues /i = -0.4,-0.6, and -0.8) and its dispersion (<t = 0.2, 0.3, 
and 0.4) for M. = 10^-^ Mrs, (Figure [Hi. We find that ^ does 
have a mild impact on the redshift position where black hole 
spin begins to decline; smaller /i results in a later decrease 
of spin. However, in comparison with Figure |7] the adopted 
uncertainty of A/i = ±0.2 does not alter the overall cosmic 
downsizing behavior for the spin evolution. Future system- 
atic studies should try to directly constrain the intrinsic Ed- 
dington ratio distribution for a given black hole mass, instead 
of inferring it from the observed Eddington ratio distribution 
combined with the AGN LF. The dispersion of the log-normal 
distribution a has only a very minor influence on the radiative 
efficiency. 

6.3.3. Compton-thiclc AGNs 

We adopt a simple, redshift-dependent correction factor to 
account for the selection bias of Compton-thick AGNs in 
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Hopkins et al.'s (2007) bolometric LFs (Equation[T3]l. Its pos- 
sible dependence on luminosity is neglected. From Equation 
it is apparent that a higher correction factor would en- 
hance the energy density and therefore lead to a larger radia - 
tive efficiency (see also iMartmez-Sansigre & Tavloiil2009l) . 
By analogy, if obscuration anti-correlates with luminosity, the 
radiative efficiency of low-mass black holes will be reduced. 
This will influence the relation between spin and mass ob- 
tained in Figure |6] Unfortunately, the overall luminosity de- 
pendence of the fraction of Compton-thick AGNs remains un- 
settled still (e.g., see Gilli et al. 2010), precluding us from ap- 
plying a rigorous treatment. Figure |9] shows the influence of 
the correction for Compton-thick AGNs on the radiative ef- 
ficiency for 1O^^M0 holes. The results are significantly af- 
fected by the normalization factor /c,o, but they are quite in- 
sensitive to the redshift parameter a. In any case, the over- 
all evolution trend of radiative efficiency is qualitatively pre- 
served. 

6.3.4. SMBH Mergers 

As in most previous works, we neglect the term for black 
hole mergers in the continuity equation and, motivated by 
Soltan's argument, assume that black hole growth is mainly 
driven by baryon accretion. As mentioned above, previ- 
ous studies show that mass accretion dominates over merg- 
ers in determini ng the mass growth and spin distributio n 
of black hole s (iBerti & Volonteril l2008t iKing et al.l l2008b . 
iShankaret al.1 (I2009L I2010h .~ using theoretically predicted 
merger rates from hierarchical structure formation models, 
also conclude that the effect of mergers is minor compared 
with mass accretion. However, based on a semi-analytic 
model of hierarchical galax y formation and evoluti on incor- 
porating black hole growth. lFanidakis et al.l (1201 Ibh recently 
showed that mergers do contribute significantly to the fi- 
nal black hole mass for high-mass black holes (e.g., M, > 
10^ Mq). Both of these calculations are subject to the uncer- 
tainties of the adopted recipes for merger rates of black holes, 
which are not well understood yet. 

6.4. Comparison with Previous Results 

Several recent parallel studies have attempted to observa- 
tionally estimate SMBH radiative efficiency or spin. Di- 
rectly starting frorn the definition of radiative efficiency, 
iDavis & Laoii (1201 Ih obtained the efficiencies of individual 
sources in a sample of 80 Palomar-Green quasars using the 
bolometric luminosity and the accretion rate determined from 
accretion disk model spectral fits. They found an average 
r/ « 0.1 and a trend of i ncreasing efficiency w ith increasing 
black hole mass (see also lRaimundo et ani201 lb . The redshift 
range of their sample is confined to z < 0.5. The correla- 
tion of the rad i ative e fficiency with black hole mass found by 
iDavis & Laod (1201 ll) is not completely compatible with our 
results (i.e. redshift bins z = 0.3 and 0.6 in Figure|6]l, although 
both calculations suffer from large uncertainties. 

Another common but highly model-dependent way to con- 
strain SMBH spin is through the jet/radio power of radio- 
loud AGNs. In a series of works, lOalyl (i2009allR IMH) 
analyzed the radio beam power using large samples of ex- 
tended radio sources and then determined the black hole spin 
with the aid of a correlation between the beam power and 
spin. Their results generally suggest that black hole spin in 
radio sources decreases moderately from redshift z w 2 to 
the local Universe. Conversely, assuming that the efficiency 



of jet production is uniquely dependen t on black hole spin, 
IMartmez-Sansigre & RawUngsl (1201 lal) inferred the spin dis- 
tribution by modeling the radio luminosity function of radio 
sources with high- and low-excitation narrow emission lines. 
They found that the best-fit spins show a bimodal distribu- 
tion, whose peaks correspond to the high- and low-excitation 
radio sources, respectively. In particular, their model fa- 
vors a trend in which the typical black hole spin increases 
slo wly toward low redshift since z ^ 3. A follow-up study 
by IMartmez-Sansigre & Rawlingsl (1201 Ibl) . based on a simi- 
lar method, gave an even steeper evolutionary trend of black 
hole spin from z « 1 to z « 0. In terms of spin evolution, our 
re sults obtained h ere ar e qualitatively consistent with those 
of lDalvl ( I2009allbl 1201 Ih . However, the stark differences be- 
tween the results of Daly and Martmez-Sansigre & Rawlings 
suggest that approaches that rely on the radio power to deduce 
black hole spin should be treated with caution. 

A number of semi-analytic models of galaxy formation 
have attempted to incorporate the growth of black holes and 
track their spin evolution. Although the prescriptions for im- 
plementing gas fueling from large to small scales are un- 
certain and differ from study to study, there is broad con- 
sensus that the final black hole spin strongly depends on 
the assumed configuration of the angular momentum of the 
accreted gas (I Volonteri et al.ll2005l: iBerti & Volonteril 120081: 
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120091: lFanidakisetal.ll20Ila|bl). Th e models of 
(120091) and lFanidakis etal.1 fcoilah predict that 
the average black hole spin increases with mass. At first 
glance, this is qualitatively consistent with our results at high 
redshifts. However, they ascribe their results to the major 
mergers of black hole s, which exclusively produ ce remnant 
holes with high spin (iHughes & Blandfordll2003h . Our cal- 
culations, by contrast, explicitly neglect black hole mergers 
and assume that mass accretion dominates black hole growth. 
We contend that this is a reasonable assumption. Indeed, 
many recent studies cast doubt on the importance of the role 
played by major mergers in the evolution of galaxies in gen- 
eral(e.g.. lNaab et al.ll2009HTavlor et alllloToriwiiEams et all 
1201 IL a nd references therein), a nd active galax ies in particular 
(e.g..lCistern as et a l. 2011 ; Ro sario et al.l20I ll and references 
therein), especially for redshifts below z ~ 2. 

7. IMPLICATIONS ON ACCRETION SCENARIOS 

Two pos sible accretion scena rios affect black hole 
spins (e.g [King & Pringle 2006 ; IBerti & Volonteril [20081: 
'King elaU 120081: IMartmez-Sansigre & Tavlod l2009h . (I) 
SMBHs grow through prolonged accretion episodes, during 
which the holes are quickly spun up to their maximum rateby 
captu ring material with constant angular momentum dThornej 
(2) SMBHs grow via many short-lived and random 
accretion episodes, acquiring gas with initially random ori- 
entations with respect to the rotation axis of the hole, which 
efficiently cancels out the net ang ular momentum and leads 
to moderately rotati ng hole s (iKing et al.l 120081 : I Wang et aP 
l2009t iLi et alllMol) . As in IWang et al.ll2009l we appeal to 
episodic, short-lived random accretion as the driving mecha- 
nism to explain the intense decline of black hole spins since 
z « 1-2. The situation may be fundamentally different at 
higher redshifts (z > 2), when galaxies are characteristically 
more gas-rich and major mergers are more prevalent. At these 
earlier epochs, it is natural to expect black hole growth to be 
more dominated by episodes of prolonged accretion, leading 
to high spins and high radiative efficiencies, particularly for 
high-mass systems. We propose that this is the fundamen- 
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tal explanation for the two regimes of spin evolution seen in 
our analysis (Figure|5]i: black holes spin up by prolonged ac- 
cretion, possibly aided by major mergers, up to z « 2, and 
thereafter spin down by episodic accretion driven by minor 
mergers and internal secular processes. 

Random accretion inevitably leaves the accretion disk mis- 
aligned with respect to the spin axis of the black hole. In this 
case, the presence of viscosity combined with Lense-Thirring 
precession will induce the inner portion of the disk to align 
or anti-align its orbital angular momentum with that of the 
black hole out to a transition radius, beyond which the disk 
retai ns its initial inclination (so -called Bardeen-Petterson ef- 
fect; lBardeen & Pettersonll975h . The characteristic extension 
of the transition radius and the time scale for alignment or 
anti-alignme nt depend on t he black hole mass and spin (e.g., 
Eringle 1992L lKinget al.l2005) . The Bardeen-Petterson effect 
may exert a vital i nfluence on the spin evolution of the black 
hole (e.g., see also lPerego et al.ll2009h . We defer a complete 
study of the relation between spin and mass during multiple 
episodes of black hole growth to a third paper of this series 
(Y.-R. Li et al. 2012, in preparation). 

8. CONCLUSIONS 

We derive the mass-dependent cosmological evolution of 
the radiative efficiency for mass accretion, which, according 
to the standard model of accretion disks, can be used as a 
surrogate to estimate the black hole spin. The calculated ra- 
diative efficiency generally increases with black hole mass, 
roughly as 77 oc Mj^'^ at high redshifts (z > 1), but the trend re- 
verses at low redshifts, such that rj increases with decreasing 
M,. High-mass black holes (M, > IO^^Mq) maintain a high 
efficiency {rj « 0.3) at z w 2, which then declines strongly to- 
ward lower redshifts. The evolutionary pattern for lower mass 
black holes is somewhat more complicated, but it is generally 
consistent with the radiative efficiency decreasing since z w 1 
in like manner Most importantly, we find that the efficiency 
of high-mass black holes begins to decline earlier than that 
of low-mass black holes, qualitatively similar to the cosmic 
downsizing of AGN activity. Assuming that the radiative effi- 
ciency provides an effective indirect measure of the black hole 
spin, we propose the following picture for the spin evolution 
of SMBHs: 



• The evolution of the spin of the black hole tracks the 
growth history of its mass and can be characterized by 
two regimes: an initial phase of mass accumulation 
from prolonged accretion that spins up the hole, fol- 
lowed by a period of random, episodic accretion that 
spins down the hole toward lower redshifts. 

• The evolution of the spin, like the global pattern of 
AGN activity, exhibits "cosmic downsizing". Relative 
to lower mass black holes, high-mass systems gain their 
masses earlier, reach the peak of their AGN activity ear- 
lier, and begin to spin down earlier. Random accretion 
dominates their evolution below z ~ 2, whereas lower 
mass holes transition to this phase later, at z ~ 1 . 

Finally, it is worth listing the principal assumptions used 
in our calculations, on which better observational constraints 
would be welcomed: (1) the scaling relation between SMBH 
mass and the mass of the bulge of the host galaxy applies 
over the redshift range < z < 2, so that the SMBH mass 
function can be inferred from the galaxy luminosity and stel- 
lar mass functions; (2) black holes gain their mass through 
baryon accretion and black hole mergers are negligible; (3) 
black hole growth mainly occurs during quasar phases char- 
acterized by universal Eddington ratio distribution; (4) the 
number of Compton-thick AGNs is comparable to that of 
Compton-thin AGNs, and they both have a similar Edding- 
ton ratio distribution as that of unobscured sources. Future 
deep multiwavelength surveys are awaited to test and refine 
these assumptions. 
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